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Abstract. - Strongly correlated two-dimensional electronic systems subject to a perpendicular 
magnetic field at lowest Landau level (LLL) filling factors: 1/2, 1/4 and 1/6 are believed to 
be composite fermion (CF) Fermi liquid phases. Even though a Bose Laughlin wave function 
cannot describe these filling factors we investigate whether such a wave function provides a 
lower energy bound to the true CF Fermi liquid energies. By using Monte Carlo simulations 
in disk geometry we compute the Bose Laughlin energies and compare them to corresponding 
results for the spin-polarized LLL CF Fermi liquid state and avalable data from literature. We 
find the unexpected result that, for filling factors v = 1/4 and 1/6, the Bose Laughlin ground 
state energy is practically identical to the true CF liquid energy while this is not the case at 
v — 1/2 where the Bose Laughlin ground state energy is sizeably lower than the energy of the 
CF Fermi liquid state. 



Modulation doped GaAs / Al x Ga\^ x As heterojunctions provide an almost ideal experimen- 
tal realization of a two-dimensional electron gas (2DEG). In these devices, the application of 
a strong perpendicular magnetic field quenches the kinetic energy and eventually brings the 
electrons in a regime where correlations are of utmost importance. In cases where the per- 
pendicular magnetic field is very strong all electrons partially occupy the lowest Landau level 
(LLL) and have their spins fully polarized in the direction of the magnetic field. The phase 
diagram of a 2DEG in a strong perpendicular magnetic field at filling factors < v < 1 is in- 
tricate with competing liquid and Wigner solid phases. At filling factors, v — 1/3, 2/5, 3/7, . . . 
and v — 1/5, 2/9, . . . electrons condense into an incompressible liquid state and the fractional 
quantum Hall effect (FQHE) occurs. [1] It is believed that at even denominator filling factors, 
v = 1/2,1/4 and 1/6 the electrons form a compressible liquid state [2-4] while for filling 
factors, < v < 1/7, Wigner crystallization occurs. [5-8] 

On the theoretical side, the accurate prediction of the critical filling factor, v c where the 
liquid-solid transition occurs is very difficult. For an ideal 2DEG system free of any disorder 
the best available theoretical estimates [5] suggest v c = 1/6.5. In experimental realm, the 
situation is more complex. Various measurements [9] on high quality samples, yet with a 
finite amount of disorder, show the existence of a reentrant solidlike phase around the FQHE 
state, v — 1/5 which likely represents an electronic solid phase pinned to fluctuations of the 
potential. The majority of theoretical studies which are based on models free of any disorder 

© EDP Sciences 



2 



EUROPHYSICS LETTERS 



agree that the critical filling factor, v c below which the Wigner crystal forms is close to and 
probably slightly larger than filling 1/7 in excellent agreement with the estimate of Lam and 
Girvin. [5]. A recent study [10] based on the exact diagonalization method found strong 
evidence that the Wigner crystal forms at filling factors, v < 1/7. Such study considered the 
filling factors, v = 1/6, 1/7, and 1/8. It was found that for v = 1/8, as well as v = 1/7 the 
translational symmetry (signature of liquid phase) is broken and the system has 2D crystalline 
order. However, the system continued to have its translational symmetry for the case of filling 
factor v = 1/6, suggesting that the value of the critical filling factor, v c is between slightly 
above 1/7, but below 1/6. 

The principal FQHE states at v = 1/3 and 1/5 are thoroughly explained and very well 
described by the Laughin wavefunction [11] while the other FQHE states at filling factors, 
v = p/(2mp + 1) ( p,m - integer) are readily understood in terms of the composite fermion 
(CF) theory. [12] The limit (p — > oo) of such FQHE states [13, 14] corresponds to even- 
denominator filled states, v = 1 /(2m) which are believed to be compressible Fermi liquid 
states qualitatively different from the FQHE states of the originating sequence. [15] The basic 
understanding now is that the low-temperature phase of fully spin-polarized electrons at filling 
fraction v — l/(2m) is a CF Fermi liquid phase. A trial wave function for the CF Fermi liquid 
ground state at v = l/(2m) has been written down by Rezayi and Read [16] and reads: 



Fermi — ±LLL 



det\e ik ^\ t> Bose , (1) 



where $B 0se is the Bose Laughlin wave function for filling factor v = l/(2m): 

*Bose = f[(zi - Zj ) 2m exp[-f;Jgf ) . (2) 



i<3 

The Rezayi-Read (RR) CF Fermi wave function is obtained after the product of the Bose 
Laughlin wave function with a Slater determinant of plane waves is fully projected into the 
LLL by means of the projection operator, Plll- In the above expressions, N is the number 
of electrons that occupy the N lowest-lying single-particle plane wave states labeled by the 
momenta {k a } consistent with an ideal 2D spin-polarized Fermi gas, zj = Xj + iyj is the 
position coordinate for the j-th electron in complex notation, lo = \ffT/ (eB) is the magnetic 
length, B is the perpendicular magnetic field and — e(e > 0) is the electron's charge. 

A straightforward motivation to study the Bose Laughlin wave function in relation to the 
true CF liquid state is justified on the naive expectation that, because of the Bose statistics, 
the Bose Laughlin energy may constitute a reasonable lower energy bound to the CF Fermi 
liquid energies for the given even-denominator filling factors. While this is an expectation 
that we cannot prove exactly, an " a-posteriori" justfication may be provided if the numerical 
results are consistent with it (although we caution that this is not a definitive proof in any 
way) . The assumption that corresponding Bose Laughlin and CF liquid ground state energies 
should be reasonably different from each other (given the very different nature of the wave 
functions) is more relevant to us in contrast to the question of whether the Bose Laughlin 
wave function is an exact lower energy bound of the true CF liquid energy. 

The main result of this work is not to assert the relevance of the Bose Laughlin wave 
function at filling factors 1/2, 1/4, and 1/6, but to report a surprising and unexpected finding, 
namely that the energy of the Bose Laughlin wave function at v = 1/4 and 1/6 (but not 1/2) 
is practically identical to the energy of the corresponding true CF Fermi liquid state. 

The model that we adopt consists of N fully spin-polarized electrons moving in a two- 
dimensional (2D) space subject to a strong perpendicular magnetic field. The electrons are 
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embedded in a uniform neutralizing background of positive charge which is spread uniformly 
within a finite disk of area fijv = i" R% and radius i?jv- The electrons can move freely all 
over the 2D space and are not constrained to stay inside the disk. The uniform density of the 
system is po = l//(2ttIq). The radius of the disk is determined from the condition: pa = N/Qn 
which gives a finite disk radius Rn — 1§\/2N jv. The kinetic energy per electron corrisponding 
to the Bose Laughlin wave function is (K)/N — hu) c /2, where uj c is the cyclotron frequency. 
Since the kinetic energy per electron is a mere constant the problem reduces to calculate the 
expectation value of the potential energy operator which consists of three terms: 

V = Vee + Vcb + Vbb , (3) 

where V ee , V e b and Vbb are, respectively, the electron-electron, electron-background and background- 
background potential energy operators. For a pure Coulomb interaction they are given by: 

N 2 N „ o „2 



(4) 

Monte Carlo (MC) simulations in disk geometry [17, 18] are used to calculate the expecta- 
tion value of the electron-electron and electron-background interaction potentials, (V ee )/N 
and (Veb) /N, respectively. The background-background interaction energy can be calculated 
exactly and is: 

(Vbb) _ 8 /WVe^ 
JV 3itV 2 ! ' [ 1 

In our MC simulations we adopt the wellknown Metropolis algorithm [19]. In this algo- 
rithm, the expectation value of any operator is estimated by averaging its value over numerous 
electronic configurations of {fi, . . . , rjv} coordinates. Our runs consist of 100,000 "equilibra- 
tion" MC steps ( a step consists of attempts to move all the electrons one by one ) and up 
to 2 x 10 6 "averaging" MC steps used to calculate the expectation value of chosen operators. 
The simulations are done for systems of JV = 4, 16, 36, 64, 100, 144, 196, 256, 324 and 400 elec- 
trons. The correlation energy per particle in the thermodynamic limit (iV — > oo) is obtained 
by fitting the finite- N energies, (V)/N to a second-order polynomial function as described in 
Ref. [17]. The resulting polynomial fits are: 

%^ = (-0.4841(5) + ^« _ g ( 

%^ = (-0.3614(03) + a03 ^ 964) - °W*t> ) g , (6) 
%^ = (-0.3013(31) + MIMHM _ om^m^ j 

The finite- N correlation energies per particle for the Bose Laughlin state at v = 1/2, 1/4 
and 1/6 and their thermodynamic limit values are displayed in Table fl] We estimated that the 
statistical uncertainty [which scales as 1 / y/MC steps] of the energy data reported in Table [T] is 
in the fifth digit after the decimal point which is rounded. Similarly, the statistical uncertainty 
of the extrapolated values is in the fifth digit after the decimal point. We also computed the 
single-particle density function, p(r) = (^2^ =1 8(r — ri)^ and the pair distribution function 

which for a liquid with uniform density po is po g(r) = (j^iLi — 1^ ~ 'j'DV ^ ne 

single-particle density of the Bose Laughlin state shows the expected nonuniformity near the 
disk boundary (see Fig. [TJ. Note that the non- uniformity near the disk edge always persists 
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Table I - Correlation energy per particle for the Bose Laughlin state at filling factors: v — 1/2, 1/4 
and 1/6. Results obtained from a standard Monte Carlo simulation in disk geometry. Energies are in 
units of e /Iq. The statistical uncertainty in the last digit of energy is shown in paranthesis. 
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Fig. 1 - One-body density function, p(r)/po, for the Bose Laughlin states v — 1/2, 1/4 and 1/6 as a 
function of the distance r/lo from the center of the disk for systems with N = 400 electrons. 

and increases as the filling factor decreases. The pair distribution function for v = 1/2, 1/4 
and 1/6 obtained after a MC simulation in disk geometry for N = 400 electrons is shown in 
FigH 

We compare the Bose Laughlin energies in the thermodynamic limit as displayed in 
last row of Table. U to respective values of the LLL-projected spin-polarized CF Fermi liq- 
uid state energies [Eq.fT])] at same filling factors. A search in the literature provides us 
with the following values for the CF Fermi liquid energy at v = 1/2: — 0.466 e 2 //o from 
Ref. [20] and —0.46557(6) e 2 /lo from Ref. [21]. Both of these values are higher than the value 
-0.4841(50) e 2 /l Q which corresponds to ^ BoS e- Ref. [20] gives the energy: -0.3608 e 2 /l for 
the v = 1/4 state, a value that is higher but very close to the energy, —0.3614(03) e 1 /Iq corre- 
sponding to the <F Bose state. We were unable to find data for the energy of the ^ Fermi state 
at v — 1/6. While there are more data for the unprojected version of the CF Fermi liquid 
wave function obtained with different methods [22-25], we found only few results for the fully 
LLL projected CF Fermi liquid states at the filling factors considered in this work. In order 
to have a larger pool of results to whom to compare the v = 1/2, 1/4 and 1/6 energies we 
pursue a different path and resort to interpolation formulas that give the energy per particle, 
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Fig. 2 - Pair distribution function of the Bose Laughlin state at filling factor v — 1/2, 1/4 and 1/6 
obtained after a MC simulation in disk geometry for N = 400 electrons. 



E(v) of the electronic liquid state as a function of the LLL filling factor: < v < 1. A crude 
empirical formula of this form has been suggested by Laughlin (L) [11] and reads: 

E L (v) = 0.814 (0.230 */>- 64 - l) \- . (7) 

Another formula proposed by Levesque, Weis, and MacDonald (LWM) [26] has the form: 

e 2 

E L wm{v) = -0.782133 y/v (l - 0.211 jA 74 + 0.012 v 1 - 7 ) — . (8) 

to 

Both expressions, Ez,{y) and E^v/Miy) are derived by fitting the energy of Laughlin wave 
functions and violate the particle- hole symmetry condition in the LLL. 

A third interpolation formula which satisfies the particle-hole symmetry condition in the 
LLL and therefore provides a more accurate interpolation for the dependence of energy on the 
LLL filling factor is given by Fano and Ortolani (FO) [27] and has the form: 



E FO {v) 



vf' 2 



■J ~v - 0.782133 V£ (1 - vf 12 + 0.55 v(l- vf - 0.463 v zj2 (1 - ,, , — 
V 8 J to 

(9) 

We use the El(i/), Elwm(v) and Epo{v) interpolation formulas to generate more data 
for the energies at v = 1/2, 1/4 and 1/6 and compare the interpolated values to Bose Laughlin 
energies at same filling factors. Such energies (in units of e 2 /lo ) are given in Table HT1 

The results in Table. [TT1 and the v = 1/2 results from Ref. [20,21] suggest that the E^(v) 
and Elwm{v) formulas are too crude to describe the energy of the electronic liquid states 
for the whole range of LLL filling factors. Clearly the Bose Laughlin energy at v = 1/2 
is distinctly smaller than the CF Fermi liquid energy of Ref. [20, 21] and FO interpolation 
energy, as well. The FO value at v = 1/2 is close to the spin-polarized CF liquid state 
value [21], although we notice that it is closer to the energy of the spin-unpolarized (singlet) 
CF state, —0.46953(7) e 2 /lo reported in Ref [21]. This is not the case for the other two filling 
factors, v = 1/4 and v = 1/6, where all the interpolation energies including the FO values are 
slightly lower than the Bose Laughlin energies, though very close to them. The more accurate 
polarized CF Fermi liquid state energy [20] (—0.3608 e 2 /l Q at v = 1/4) is larger than the Bose 
Laughlin energy, though it is important to note that the two values are so close that they can 
be considered practically identical. 



6 



EUROPHYSICS LETTERS 



Table II - Energy per particle for the LLL states v — 1/2, 1/4 and 1/6 obtained from interpolation 
formulas: Eq. (\7ty, Eq. and Eq. as compared with Bose Laughlin state energies. Energies are 
in units of e 2 /lo- The statistical uncertainty in the last digits of energy is shown in paranthesis. 
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The fact that interpolation energies are slightly lower than the Bose Laughlin energies 
at v = 1/4 and 1/6 surely does not invalidate the possibility of the Bose Laughlin wave 
function being a lower energy bound to the true CF Fermi liquid energy since it is well 
known [28] that the true energy at any filling factor (excluding the Laughlin filling factors) 
is always higher than the approximate energy obtained from L, LWM and FO interpolation 
formulas that are smooth functions of v and should be handled with care. For electronic 
LLL Laughlin wave functions (Fermi or Bose) and a bare Coulomb potential, the Haldane 
pseudopotential parameters, V m are known to be all positive and decrease monotonically as 
the relative angular momentum, m increases. In contrast to the electronic case, the Haldane 
pseudopotential parameters (for both bare and screened Coulomb potentials) corresponding 
to the CF Fermi liquid states in Eq. ( [1]) show a distinct non-monotonic behavior with some 
of them even becoming negative for large values of relative angular momentum. [29] In the 
CF language, a state with filling factor, v = p/(2mp+ 1) represents a CF system with p filled 
CF LL-s which weakly interact with each other by some residual interaction. In such a case, 
the CF Haldane pseudopotentials (both small p and large p) show a qualitatively different 
behavior from the electronic LLL Laughlin pseudopotentials. Such a qualitatively different 
behavior is likely associated with the weakly interacting nature of the CF-s as compared to 
pure electrons. [30] Because of the different nature of the two liquid wave functions, [Eq.([T]) 
and Eq.([2])] and their respective Haldane pseudopotentials, it was unexpected to find that 
the Bose Laughlin and the true CF Fermi liquid energies are practically identical for filling 
factor v = 1/4 and 1/6, while clearly distinct for v = 1/2. 

In summary, we used MC simulations to calculate the ground state properties of Bose 
Laughlin wavefunctions at compressible filling fractions v — 1/2, 1/4 and 1/6. The calculated 
ground state energies are compared with results for spin-polarized CF wavefunctions for these 
filling factors. We found that the Bose Laughlin wavefunction ground state energy is lower 
than the CF Fermi liquid energy for v = 1/2, whereas for other filling factors, v = 1/4 and 
1/6 we find the surprising result that the Bose Laughlin ground state energy is practically 
identical to the energy of the true CF liquid state as calculated from different studies and 
from estimates obtained from various interpolation schemes. Because Bose Laughlin and true 
CF liquid energies are very different at v = 1/2 we naively expected to see a similar tendency 
apply to the other two even-denominator- filled CF liquid states, v = 1/4 and 1/6, however 
unexpectedly we found out that the Bose Laughlin and true CF liquid energies are practically 
identical at filling factors 1/4 and 1/6. The physical significance and robustness of this finding 
is not yet clear to us and further work is needed to address this issue. 
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